Spin relaxation in n-type GaAs quantum wells from a fully microscopic approach 
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We perform a full microscopic investigation on the spin relaxation in n-type (001) GaAs quantum 
wells with Alo.4Gao.6As barrier due to the D'yakonov-perel' mechanism from nearly 20 K to the 
room temperature by constructing and numerically solving the kinetic spin Bloch equations. We 
consider all the relevant scattering such as the electron-acoustic-phonon, the electron-longitudinal- 
optical-phonon, the electron-nonmagnetic-impurity and the electron-electron Coulomb scattering 
to the spin relaxation. The spin relaxation times calculated from our theory with a fitting spin 
splitting parameter are in good agreement with the experimental data by Ohno et al. [Physica E 6, 
817 (2000)] over the whole temperature regime (from 20 K to 300 K). The value of the fitted spin 
splitting parameter agrees with many experiments and theoretical calculations. We further show 
the temperature dependence of the spin relaxation time under various conditions such as electron 
density, impurity density and well width. We predict a peak solely due to the Coulomb scattering 
in the spin relaxation time at low temperature (< 50 K) in samples with low electron density {e.g., 
density less than 1 x 10^^ cm~^) but high mobility. This peak disappears in samples with high 
electron density {e.g. 2 x 10^^ cm^^) and/or low mobility. The hot-electron spin kinetics at low 
temperature is also addressed with many features quite different from the high temperature case 
predicted. 



PACS numbers: 72.25.Rb, 72.20.Ht, 71.10.-w, 67.57.Lm 



I. INTRODUCTION 



Much attention has been devoted to the elec- 
tron spin dynamics in semiconductors for the past 
three decades. "'^'^ Especially, recent experiments have 
shown extremely long spin lifetime (up to hundreds of 
nanoseconds) in n-type bulk Zinc-blende semiconduc- 
tors (such as GaAs).'^''*'^ Moreover, a lot more investi- 
gations have been performed on various low dimensional 
Systems,6.7.8,9,i0,ii,i2,i3,i4.i5,i6.i7,i8.i9,20 spin lifetime 

as long as tens of nanoseconds has been reported in 
(llO)-oriented GaAs quantum wells (QWs)"'^'''^^ at room 
temperature and inp-type GaAs:Be/Al2;Gai_a;As double 
hetero-structures^ at low temperature. In these stud- 
ies, understanding the spin relaxation/dephasing (R/D) 
mechanism is one of the most important problems as it is 
the prerequisite for the application of the spintronic de- 
vices. It is understood that the D'ayakonov-Perel' (DP) 
mechanism is the leading spin R/D mechanism in n-type 
Zinc-blende semiconductors.^^ This mechanism is com- 
posed of the contribution from the Dresselhaus term,^^ 
which is due to the lack of inversion symmetry in the 
Zinc-blende crystal Brillouin zone (sometimes referred to 
as the bulk inversion asymmetry), and that from the 
Rashba term,^^ which originates from the asymmetric 
potential within a QW along the growth direction (some- 
times referred to as the structure inversion asymmetry). 
Both appear as effective magnetic fields. For narrow 
(001) GaAs QW without the additional large bias volt- 
age, the Dresselhaus term is the leading term:^'*'^^ 
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in which (fc^) represents the average of the operator 
— {d/dzY over the electronic state of the lowest subband. 
Under the finite square well assumption. 



{kl) 



AA, 
— 2/3- 



cos^ ^ 



) , 



(4) 



where f and (3 are the lowest energy solutions of the equa- 
tions 
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^ with V , a and m* denoting the 
well depth, well width and the effective mass respectively. 
In the limiting case, \\mv—>oo{k1) = {ir/a)'^. 7 is the spin 
splitting parameter.^ There are a lot of theoretical in- 
vestigations on the spin R/D due to the DP mechanism 
lately.2^'2^'2^ Most of them are within the framework of 
single-particle approximation and the Coulomb scatter- 
ing is thought to be irrelevant in the spin R/D. 

Recently Wu et al. performed a full microscopic 
investigation on the spin R/D and showed that the 
single-particle approach is inadequate in accounting for 

the spin R/D.29.30,31,32,33,34,35,36,37,38,39,40 ^i^jg ^p- 

proach, the momentum dependence of the effective mag- 
netic field (the DP term) and the momentum depen- 
dence of the spin diffusion rate in the direction of 
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the spacial gradient^* or even the random spin-orbit 
interaction^^ serve as inhomogeneous broadening. ■^"■'^^ 
In the presence of the inhomogeneous broadening, any 
scattering (even the spin-conserving scattering), includ- 
ing the Coulomb scattering, ^^'^^'^^'^^ can cause irre- 
versible dephasing. Moreover, this approach also in- 
cludes the counter effect of the scattering to the inho- 
mogeneous broadening, i.e., the suppression of the inho- 
mogeneous broadening by the scattering. Finally, this 
approach is valid not only near the equilibrium, but 
also far away from the equilibrium, ^^'^^ and is appli- 
cable to both the strong (|r2|Tp <^ 1) and the weak 
([r2|Tp » 1) scattering hmits,^^'^" with Tp representing 
the momentum relaxation time. In the weak scatter- 
ing limit, the counter effect of the scattering is less im- 
portant and adding additional scattering (including the 
Coulomb scattering) causes stronger spin R/D. Whereas 
in the strong scattering limit, adding additional scat- 
tering always increases the spin R/D time. The fea- 
ture is more complicated when \fl\Tp ~ 1.'^^ In above 

Studies,29'30,31,32,33,34,35,36,37,38,39.40 j^^ve been foCUS- 

ing on the high temperature regime (T > 120 K) where 

the electron-acoustic (AC) phonon scattering, which 
is more complicated in numerical calculation than the 
electron-longitudinal optical (LO) phonon scattering, is 
negligible. In this paper, we extend the scope of our 
approach to study the spin kinetics at low temperature 
regime by including the electron- AC phonon scattering. 
Moreover, we compare the spin relaxation time (SRT) 
obtained from our theory with the experimental data 
over a wide temperature regime and show the excellent 
agreement of our theory with the experiment. We fur- 
ther show that the Coulomb scattering is important to 
the spin R/D not only at high temperatures, ^^'^s. 39 j-^^^ 
also at low temperatures. The electron density, impu- 
rity density, well width, temperature and electric field 
dependences of the SRT are studied in detail. 

This paper is organized as follows: In Sec. II we set 
up the model and give the kinetic spin Bloch equations. 
In Sec. Ill we compare our results with the experimen- 
tal data. Then, we investigate the temperature depen- 
dence of the spin relaxation under different conditions 
such as electron densities, impurity densities and well 
widths in Sec. IV. The effect of Coulomb scattering is 
also addressed. The hot-electron effect in spin relaxation 
is investigated in Sec. V. We summarize in Sec. VI. 



II. KINETIC SPIN BLOCH EQUATIONS 

We start our investigation from an n-type GaAs (001) 
QW with the growth direction along the ^-axis. A mod- 
erate magnetic field B is applied along the x-axis (in the 
Voigt configuration). The kinetic spin Bloch equations 
can be constructed by using the nonequilibrium Green 
function method:*^ 

Pk,aa' — eE • VkPk,<Tt7' = Pk,<7<7' \coh + Pk,<7<7' \scatt^ (6) 



with pk.crcr' representing the single particle density ma- 
trix elements. The diagonal and off-diagonal elements 
give the electron distribution functions /kcr and the spin 
coherence Pk,a-cr- The second term in Eq. (6) de- 
scribes the energy input from the external electric field 
E. The coherent terms pk,crcr'\coh describe the pre- 
cession of the electron spin due to the applied mag- 
netic field B and the effective magnetic field Jl(k) 
[Eqs. (1-3)] as well as the effective magnetic field from 
the Hartree-Fock Coulomb interaction. pk,aa'\scatt in 
Eq. (6) denote the electron-LO-phonon, the electron- 
AC-phonon, the electron-nonmagnetic impurity and the 
electron-electron Coulomb scattering. Their expressions 
are given in detail in Rcf. 36, except the additional 
matrix elements of the electron- AC-phonon scattering. 
For the electron-AC-phonon scattering due to the de- 
formation potential, the matrix elements are given by 

■9q de/ ^ 2dv'^i [-^(^g-s ) P ^-^d for the scattering due 
to the piezoelectric coupling, the matrix elements read 
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\I{iqz)\^ for the longitudi- 



ah, 

nal phonon and 5^ ^^ = l^J^i^t^'v + ^U't + 



(3q.qyg.) for the transverse phonon.*^ 

Here Q = (q, q^); Q = ^ql+ql + qf, S = 8.5 eV is the 

deformation potential; d = 5.31 g/cm^ is the mass den- 
sity of the crystal; Vsi = 5.29 x 10^ m/s {vst = 2.48 x 10^ 
m/s) is the velocity of the longitudinal (transverse) sound 
wave; kq ~ 12.9 denotes the static dielectric constant; 
and ei4 = 1.41 x 10^ V/m represents the piezoelectric 
constant. The AC phonon spectra uq\ are given by 
^Qi = VsiQ for the longitudinal mode and LUQt = VstQ 
for the transverse mode. The form factor is 
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with y = qza/2. The numerical schemes of the electron- 
electron Coulomb, the electron- impurity as well as the 
electron-LO phonon scattering have been given in detail 
in Ref. 36, whereas the numerical scheme for the electron- 
AC-phonon scattering is presented in Appendix A. The 
electron-interface-phonon scattering is negligible due to 
the thick GaAlAs barrier. In addition, as we are going 
to explore the spin R/D over a wide range of electron 
densities, in the present paper we use the screening under 
the random phase approximation^^ rather than the one 
in the limiting (degenerate or nondegenerate) cases for 
the screened Coulomb potential:^'^ 



v„ = 
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where vq = 47re^/(5^ is the bare Coulomb potential and 



P(l)(q) = /k+qg /kg- 
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In this way, we also take into account the hot-electron 
effect on the screening. 

By numerically solving the kinetic spin Bloch equa- 
tions with all these scattering explicitly included, one is 
able to obtain the spin dephasing and relaxation times 
from the temporal evolutions of the spin coherence pk^a-a 
and the electron distribution functions /k.o-. The ir- 
reversible spin dephasing time can be obtained by the 
slope of the envelope of the incoherently summed spin 
coherence^^ p = |Pk,Ti(OI> ^'^d the SRT can be de- 
fined by the slope of the envelope of the difference be- 
tween n-\ and nj^, with = /k,o- 




T(K) 

FIG. 1: (a) SRT r vs. temperature T for GaAs QW with 
a = 7.5 nm and electron density n = 4 x lO'^" cm~^ at three 
different spin-splitting parameters. Dots: experiment data; 
Dot-dashed curve: 7 = O.970; Solid curve: 7 = 70; Dashed 
curve: 7 = l.l7o- (b) Hall mobility /iHaii vs. temperature T 
(Ref. 12). 



III. COMPARISON WITH EXPERIMENT 

First, we compare the SRT obtained from our micro- 
scopic approach with the experimental data by Ohno et 
al. in 60 periods of GaAs QWs separated by 10 — 12 nm 
thick Alo.4Gao.6As barriers. The well width of each QW 
a ~ 7.5 nm and the electron density n = 4 x 10^° cm"^.-'^^ 
The well depth of the GaAs well confined by Ala,Gai_a;As 
is roughly estimated to be 65 % of 1.087a; + 0.438a;2 eV, 
for X = 0.4 and is therefore Vo = 328 mcV."*^ Differently 
from our previous fit'^^ with the experiment data at high 



temperatures by Malinowski et al.^^ where we had two 
fitting parameters, i.e., the spin splitting parameter 7 
and the impurity density Ui due to the absence of mo- 
bility data, here we have only one fitting parameter 7 
and the temperature sweeps from the very low temper- 
ature to the room temperature. The corresponding Hall 
mobilities /inaii in the experiments^ can be found in Ref. 
27, also plotted in Fig. 1(b). From the Hall mobility, one 
can deduce the impurity density by calculating the trans- 
port mobilities''^ utr = MHaii/j^Haii with rnaii = 1 for the 
electron- AC phonon scattering due to the deformation 
potential; rnaii = 7/5 for the electron- AC phonon scatter- 
ing due to the piezoelectric coupling and the electron-LO 
phonon scattering; and rnaii = 1 for the electron-ionized 
impurity scattering. "''^ 

The only fitting parameter 7 is around 



7: 



(4/3)(m7me.)(l/^2m*3£^)(,,/7rr^) , (10) 



in which 77 = A/ {Eg + A); Eg denotes the band gap; A 
represents the spin-orbit splitting of the valence band; 
m* stands for the electron mass in GaAs; and nicv is 
a constant close in magnitude to the free electron mass 
mo.^^ For GaAs when m.cv = mo, 7 = 70 = 11-4 eV-A^. 
The initial spin polarization P is assumed to be 2.5 % for 
weak polarization throughout the paper. In Fig. 1(a), the 
SRTs T obtained from our approach are plotted against 
the temperature with all the scattering included. B and 
E are taken to be zero, as in the experiment. 7 = 0.9, 
1.0 and I.I70 correspond to nicv = 1-1, 1-0 and 0.91mo 
with nicv being the only not-fully-determined parame- 
ter in Eq. (10). One finds good agreement between our 
theory and the experiment data almost over the whole 
temperature regime. When T is below 13 K, there is no 
theoretical data due to the lack of experimental data for 
the Hall mobility. Kainz et al. also fitted the same exper- 
iment data by using the single-particle theory without the 
Coulomb scattering.^'' They used a fourteen-band model 
to calculate the spin-orbit coupling. Unlike our theory, 
their results can only give the boundary values of the SRT 
in several cases rather than the exact data.^^ This is be- 
cause they did not take the full microscopic calculation, 
and the single-particle theory is inadequate in accounting 
for the spin R/D. 

The best-fitted value of the spin-orbit coupling pa- 
rameter, i.e., 7o = 11.4 eV-A^, is close to the value 
calculated by Kainz et al. {'^ 16.5 eV-A"^) using the 
multiband envelope-function approximation.^^ It is noted 
that the value for 7 in GaAs is still in debate. Usu- 
ally reported experimental values for 7 (25 — 30 eV-A^) 
in bulk material are deduced from the DP spin relax- 
ation mechanism within the framework of the single- 
particle approximation, where the Coulomb scattering is 
not included. Furthermore, the Raman scattering ex- 
periment showed that 7 = 16.5 ± 3 eV-A^ in asymmet- 
ric QW;^^ and the Hanle effect experiment showed that 
7 = 12.6 eV-A*^."''"' Theoretically, semi-empirical param- 
eterized 16 X 16 k-p calculations show that 7 = 14.9 
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eV-A'^;^^ and the self-consistent ab initio calculations 
predict 6.4 and 8.5 eV-A''.'''*' Our fitting result supports 
the last four experimental and theoretical results. It is 
also noted that our result further confirms the analytical 
result Eq. (10) obtained from the perturbation,^^ with 
nicv = mo- 



IV. TEMPERATURE DEPENDENCE OF SRT 

We now study the temperature dependence of the spin 
relaxation in detail. In the calculation, the electric field 
£^ = 0, the magnetic field B = T and the spin splitting 
parameter 7 = 70. 

We plot in Fig. 2 the temperature dependence of the 
SRT of GaAs/Alo.4Gao.6As QWs with a = 7.5 nm at dif- 
ferent impurity densities when the electron densities are 
low (n = 4x10^° cm-2) [Fig. 2(a)], medium (n = IxlO" 
cm-2) [Fig. 2(b)] and high (n 2 x lO" cm'^) [Fig. 
2(c)] respectively as solid curves. For the well width 
and the electron densities here, the linear terms in the 
DP terms [Eqs. (1-3)] are dominant, and only the low- 
est subband is relevant when T < 300 K. It is seen 
from the figure that adding impurities always increases 
the SRT. This is understood that the criterion of strong 
scattering \fl\Tp <C 1 is satisfied here at all tempera- 
tures, and therefore adding additional scattering always 
increases the SRT.^^ It is noted that Tp here has been 
extended to include r^^, i.e., the contribution from the 
Coulomb scattering."'''^ It is interesting to note; that un- 
like our previous works focusing on high temperatures 
(T > 120 K), 32.33,34,35,36,37,38,39,40 situation is more 

complicated at low temperatures. At low/medium elec- 
tron densities [Fig. 2(a)/(b)], the SRT presents a peak 
at very low temperature (near 20~30 K)/low tempera- 
ture (around 41 K) and a valley around 120 K; whereas 
at high electron density [Fig. 2(c)], the SRT increases 
monotonically with T. 

It is noted that at very low temperature (around 20 K) 
the electron-AC phonon scattering is negligible. t^'^ 
from the electron-AC phonon scattering is around 25 
ps, two orders of magnitude larger than t^'^ from the 
electron-electron Coulomb scattering. In addition, 
from the impurity scattering is around 2 ps, one or- 
der of magnitude larger than Tp^, and has a very weak 
temperature dependence. Therefore, the appearance 
of the peaks in Fig. 2(a) originates from the electron- 
electron Coulomb scattering which dominates the scat- 
tering process. Moreover, Tp^ is a nonmonotonic func- 
tion of temperature: r^^ oc T^^ at low temperature (de- 
generate limit) and r^'' cx: T at high temperature (non- 
degenerate limit). The minimum of Tp^ corresponds 
to the crossover from the degenerate limit to the non- 
degenerate one at ~ Ep/kB- T'c ~ 17 K when 
n = 4 X 10^° cm~2^ in good agreement with the peaks 
obtained from our calculation with the exact Coulomb 
scattering. Therefore the SRT increases/decreases with 



10" 



10' 







(a) 






ii='i X lu /cm 















so 100 ISO 200 2S0 300 
T(K) 



10" 



10' 




so 100 ISO 200 2S0 300 
T(K) 




SO 100 ISO 200 2S0 300 
T(K) 

FIG. 2: SRT r vs the temperature T with well width a = 7.5 
nm and electron density n being (a) 4 x 10^° cm~^, (b) 1 x lO^'^ 
cm~^ and (c) 2 x 10^^ cm~^ respectively. Solid curve with 
triangles: Ui = n; Solid curve with dots: rii = O.ln; Solid 
curve with circles: rn = 0; Dashed curve with dots: = O.ln 
and no Coulomb scattering. 



the temperature (and the Coulomb scattering) in the de- 
generate/nondegenerate regime. Once T > 120 K, the 
electron-LO phonon scattering becomes comparable with 
the Coulomb scattering and strengthens so rapidly with 
temperature that it completely surpasses the weak tem- 
perature dependence of the Coulomb scattering: t^^ 
from the electron-LO phonon scattering varies from sev- 
eral picoseconds at 120 K to several tenths of picosecond 
at 300 K, and r^*^ varies from 1 ps to several picoseconds. 
Therefore the SRT increases with T. When the electron 
density is 1 x 10^^ cm^^, Tc is nearly 41 K. Around this 
temperature, the electron-AC phonon scattering cannot 
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be overlooked although is still roughly one order of 
magnitude larger than t'^'^ . Therefore, the reduction of 
the Coulomb scattering after Tc can be partly compen- 
sated by the increase of the electron-phonon scattering. 
As a result, one can see that the decrease of the SRT af- 
ter Tc in Fig. 2(b) is much slower than that in Fig. 2(a). 
However, when the electron density is high enough, say 
2 X 10^^ cm~^ in Fig. 2(c), Tc is nearly 83 K, much larger 
than the case of low electron density. At this tempera- 
ture, the electron-phonon scattering becomes comparable 
to the Coulomb scattering and the strengthening rate 
of phonon scattering around this temperature is large 
enough to completely compensate, and even surpass, 
the weakening rate of the Coulomb scattering. Conse- 
quently the total scattering increases monotonously with 
T. Therefore the SRT increases monotonically with T. 

We further show the effect of the Coulomb scattering 

on the spin relaxation. This was first proposed by Wu 
and Ning based on the inhomogeneous broadening in- 
duced by the energy dependence of the g-factor.'^" Then, 
we used our full microscopic approach and showed that 
the Coulomb scattering makes marked contribution to 
the spin R/D when T > 120 K when the inhomogeneous 
broadening is induced by the DP term.'^^'^^''^^ At low tem- 
perature (T < 120 K), Glazov and Ivchenko used pertur- 
bation method to show that the second-order Coulomb 
scattering causes the SRT.^'' In the perturbation ap- 
proach, the Coulomb scattering contributes marginally to 
the spin R/D at high temperature. In our calculation, we 
include the Coulomb scattering to all orders of the bub- 
ble diagrams as well as the counter effect of the Coulomb 
scattering to the inhomogeneous broadening. In Fig. 2 
by plotting the SRT for the case of rii = O.ln, but with- 
out the Coulomb scattering, as dashed curves, we show 
that the Coulomb scattering makes marked contribution 
to the spin R/D over the whole temperature regime by 
increasing the spin R/D time.^'^ It is further seen from 
Figs. 2(a) and (b) that the peak disappears without the 
Coulomb scattering. This is consistent with the previous 
discussion. 

It is interesting to see that in the absence of the 
Coulomb scattering, the criterion for strong scattering 
regime |n|Tp ^ 1 is satisfied only when T > 120 K. 
Therefore the SRT increases with T when T > 120 K. 
When T < 120 K, \ fl\Tp is slightly smaller than 1, which 
is the intermediate scattering regime. The variation of 
the SRT depends on the competition between the in- 
crease of the inhomogeneous broadening and the increase 
of the scattering with the temperature.^^ For low/high 
electron density case, the temperature dependence of the 
electron- AC phonon scattering is less/more effective and 
the SRT decreases/increases with T. 

It is noted that in order to see the peaks at low electron 

density, it is important to have a high mobility sample 
(low impurity density) . This is because the ascendancy of 
the Coulomb scattering can be impaired when the impu- 
rity scattering gets large enough and the total scattering 
is mainly determined by the impurity scattering. As the 
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FIG. 3: SRT r vs. temprature T at a = 7.5 nm (curves 
with dots) and 15 nm (curves with triangles). Solid curves: 
Hi = O.ln; Dashed curves: rii = n. (a) n = 4 x 10^" cm~^ 
and (b) n = 2 X lO" cm-^. 

electron-impurity scattering depends weakly on the tem- 
perature, the temperature dependence of the inhomoge- 
neous broadening from the DP term becomes the only 
variable element. Therefore the SRT decreases monoton- 
ically with T as the solid curves with triangles in Fig. 
2 for the case of = n. This condition is not satisfied 
in the experiment by Ohno et al.,^^ and this is the rea- 
son why there is no peak in Fig. 1. However, apart from 
the peaks which are not observed yet, both the experi- 
ment and calculation show that the SRT decreases with 
temperature at low electron densities when T < 120 K. 
The SRT at high electron density increases monotoni- 
cally with temperature when the impurity density is low, 
which is also in agreement with the latest experiment by 
Harley et al.^^ 

Finally we investigate the well width dependence of the 
SRT. In Fig. 3 we plot the SRT versus temperature at 
well widths a = 7.5 nm (solid curves) and 15 nm (dashed 
curves) respectively. We choose low and high impurity 
densities n, = O.ln (curves with dots) and rii = n (ciuves 
with triangles) as well as low and high electron densities 
n = 4 X 10^° cm-2 (a) and n = 2 x 10^^ cm^^ j^^). It 
is noted that the SRT is enhanced by increasing the well 
width as {k1) in the DP term decreases with the increase 
of a. Moreover, as impurities further enhance the SRT, it 
reaches to several nanoseconds at very low temperatures 
at high impurity density. 
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V. ELECTRIC FIELD DEPENDENCE OF SRT 

We now turn to investigating tlie hot-electron effect 
on spin relaxation at low temperature. An electric field 
is applied parallel to the QW. Similar to our previous 
investigation,^^ electrons obtain a center-of-mass drift ve- 
locity (and consequently an effective magnetic field pro- 
portional to the electric field) and are heated to a tem- 
perature Te higher than T. The numerical schemes of 
solving the hot electron problem has been laid out in de- 
tail in Ref. 36.*''' We plot the electric field dependence of 
the SRT with a = 7.5 nm and T = 50 K for different im- 
purity densities at low/high electron density [n = AxlQi^^ 
cm~2/n = 2 X 10^^ cm"^) in Fig. 4(a)/(b). In the calcu- 
lation, the magnetic field _B = 4 T and the spin splitting 
parameter 7 = 7q as in the previous section. It is seen 
from the figure that unlike the high temperature case in- 
vestigated before''^ (and also see Fig. 4(c) for T = 120 K) 
where the electric field can be applied easily to around 
1 kV/cm, at low temperatures it can be applied only to 
a very small value due to the "runaway" effect. This 
is because at low temperature, the efficient electron-LO 
phonon scattering is missing and electrons are therefore 
very easily driven to very high momentum states by a 
very small electric field. 

It is interesting to see from the figure that differing 
from the high temperature case where the SRT increases 
with the electric field (see Fig. 4 (c) and also Ref. 36), 
here, for the case of low electron densities, the SRT de- 
creases with the field and for the case of high electron 
densities, the SRT decreases/increases with the field at 
high/low impurity densities. 

These features at low temperature T can be under- 
stood from the joint effects of the electric field E to the 
scattering strength and the inhomogeneous broadening 
due to the DP term. On one hand, when the electric field 
is small, the ionized-impurity scattering, whose strength 
decreases slightly with the electron temperature Tg, is 
dominant. When the electric field is further increased, 
Te and therefore the electron-AC phonon scattering is 
raised. If the impurity scattering is not too high, the 
electron-AC phonon scattering can then be dominant as 
discussed decades ago in Ref. 65. These can be seen 
from the mobilities /U = ^k/k(7/(rn*n£') obtained 
from our calculation, which are plotted in the same fig- 
ure for all the corresponding cases. One can see that 
fj, increases slightly and monotonically with E for the 
impurity-scattering-dominant case such as Ui = n; It de- 
creases monotonically with E for the case of very low/no 
impurity scattering such as rij = 0, as the electron-AC 
phonon scattering always increases with the electron tem- 
perature Te, for the case of low impurity scattering such 
as rii = O.ln, /z first increases slightly then decreases 
with E which shows the transition from the impurity- 
scattering-dominant regime to the electron- AC-phonon- 
scattering-dominant regime^^ (unless the runaway effect 
blocks the system to the later regime as shown in Fig. 
4(a) for the case of low electron densities). On the other 
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FIG. 4: SRT T (solid curves) and mobility /i (dashed curves) 
vs. electric field E. (a) T = 50 K and n = 4 x 10^° cm"^; 
(b) T = 50 K and n = 2 X ll" cm-^; (c) T = 120 K, 
n — 4x10^" cm~^ and2xll^^ cm~^ respectively. Curves with 
open circles: = n; with dots: Ui = O.ln; with triangles: 
rii = 0. Note the scales of the mobility fi are on the right side 
of the figures. 



hand, electrons are driven to the higher momentum states 
by the electric field and experience a larger effective mag- 
netic field [Eqs. (1-3)]. Therefore the inhomogeneous 
broadening is increased. This tends to decrease the SRT. 
In order to show the electric field dependence of the in- 
homogeneous broadening, we plot in Fig. 5 the electron 
temperature Tg as a function of electric field E for all the 
corresponding cases in Figs. 4(a) and (b). It is obvious 
from Fig. 5 that the increase of Tg for the low electron 
density case is much faster than that for the high one. 
Therefore, the increase of the inhomogeneous broadening 
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FIG. 5: Hot electron temperature Te vs. electric field E when 
T = 50 K. solid curves: n = 4 x 10^" cm~^; dashed curves: 
n = 2 X 11^^ cm~^. Curves with open circles: = n; with 
dots: rii = O.ln; with triangle: rii = 0. 

is the leading contribution in comparison with the elec- 
tric field effect on the scattering. Consequently, the SRT 
decreases with E for the case of low electron densities as 
shown in Fig. 4(a). For the case of high electron densities, 
when the impurity density is high such as rii = n in Fig. 
4(b), both the slight decrease of the scattering and the 
increase of the inhomogeneous broadening tend to sup- 
press the SRT. When the impurity density is low/zero, 
the strengthening of the scattering is dominant (as shown 
in Fig. 4(b) the decrease of mobility with E) in compar- 
ison with the increase of the inhomogeneous broadening. 
This makes the SRT dcicrcasc with E. Finally we point 
out that the Coulomb scattering plays an essential role 
in the spin R/D in the presence of the electric field. It 
determines the hot-electron temperature which con- 
trols the inhomogeneous broadening and the scattering 
strengths. Moreover, the Coulomb scattering itself also 
contributes to the spin R/D. 

For comparison, we also plot the SRT versus electric 
field at high temperature T = 120 K for both low and 
high electron densities with rii = O.ln in Fig. 4(c). At 
this temperature, the electron-LO phonon scattering is 
dominant. Therefore always decrease with E. For the 
well width and electron density we study, the linear DP 
term is dominant and the increase of the scattering is 
more important. Therefore the SRT increases with E. 
For QWs with larger well width so that the cubic term 
is dominant, the SRT can decrease with E as reported in 
our previous work at high temperatures.^'' 



VI. SUMMARY 

In summary, we have investigated the temperature de- 
pendence of the SRT for n-typc GaAs (001) QWs with 
small well widths from a full microscopic approach by 
constructing and numerically solving the kinetic spin 
Bloch equations with all the relevant scattering explic- 
itly included. In contrast to our previous studies at high 



temperatures (T > 120 K), we include the electron- AC 
phonon scattering which is absent in our previous stud- 
ies so that we may extend the scope of our approach to 
the low temperature regime (T < 120 K). Good agree- 
ment with experiment data^^ is obtained from our theory 
over almost the whole temperature regime by using only 
one fitting parameter 7 whose value agrees with many 
experimental and theoretical results. We show that the 
Coulomb scattering plays an essential role in spin R/D 
over all the temperature regime. 

For QWs with low electron densities but high mobil- 
ity (i.e., low impurity density), the spin R/D is mainly 
controlled by the electron-electron Coulomb scattering 
when T < 70 K. We predict a peak in the t-T curve. 
The closer the peak approaches the high temperature 
limit, the smoother the peak appears. After the peak, 
the SRT increases with temperature. Finally, such a 
peak disappears at sufficient high electron density where 
the SRT increases monotonically with temperature. We 
point out that the peak origins from the Coulomb scatter- 
ing. Specifically, it origins from the different temperature 
dependences of the Coulomb scattering at the degenerate 
and the nondegenerate limits with the transition temper- 
ature Tc ~ Ep/ku- For low electron densities, Tc < 30 K 
where the electron-phonon scattering is negligible. Then 
one may observe an abrupt peak around Tc. For medium 
electron densities, 30 K< < 70 K where the increase of 
the electron- AC phonon scattering partially compensates 
the decrease of the Coulomb scattering when T increases, 
one may observe a smooth peak around T^. Nevertheless, 
for high electron densities, Tc > 70 K, the increase of the 
electron-phonon scattering completely compensates the 
decrease of the Coulomb scattering when T rises. Con- 
sequently the peak disappears. 

At high temperature (T > 120 K) and low impurity 
density, when the well width is small so that the cu- 
bic terms in the DP terms are unimportant, the increase 
of electron-LO phonon scattering surpasses the increase 
of inhomogeneous broadening with temperature, so that 
the SRT increases with temperature. However, when the 
impurity density is so high that electron-impurity scat- 
tering is the dominant scattering mechanism, the SRT 
decreases monotonically with temperature for any elec- 
tron density. This is because the temperature depen- 
dence of the electron-impurity scattering is very weak 
and the increase of the inhomogeneous broadening with 
temperature dominates the temperature dependence of 
the SRT. We also show that larger well width leads to a 
slower spin relaxation. Moreover, in the strong scatter- 
ing limit, higher impurity density also leads to a slower 
spin relaxation. Both effects can make the SRT as long 
as nanoseconds at very low temperatures. 

The effect of electric field (i.e., the hot electron effect) 

on the spin relaxation is also investigated. We show that 
the electric field dependence of the SRT at low temper- 
ature appears again quite differently from that at high 
temperature due to the absence of electron-LO phonon 
scattering. Moreover, we further show different electric 
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field dependences of tlic SRT at low and high electron 
densities. At low electron densities, the SRT decreases 
with the electric field. When the electron density is high, 
it decreases/increases with the electric field for the case 
of high/low impurity densities. These features are in gi- 
ant difference from the high temperature case where the 
SRT increases monotonically with electric field for the 
same QWs. More experiments are needed to explore the 
predictions presented in this manuscript. 
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APPENDIX A: NUMERICAL SCHEME FOR 
ELECTRON-AC PHONON SCATTERING 



The electron-AG phonon scattering terms can be 
rewritten as 



k,(T 



dt 



dt 



AC 



AC 



{-271- 5q9,a'5(ek - ek-q - ^qq,x)[Niek - ek-q)(/k,<T - h-q,a) + /k,a(l " /k-q,, 

-Re(pkpj;_q)]}-{k^k-q} , 

{'^ m ^qg.A'^lek - £k-q " ^^qg, a) [Pk-q(/k,T + /k,i) + (/k-q,T + /k-q,i " 2)pk 

-2N{ek - ek-q)(pk - Pk-q)]| - |k <m. k - q| , 



(Al) 



(A2) 



with pk = Pk,Ti and {k k — q} standing for the 
same terms as the previous {} but with the interchange 
k ^ k-q. iV(ek-ek-q) = [exp{/3(ek-ek-q)}-l]~^ rep- 
resents the Bose distribution. The division of the trun- 
cated two-dimensional momentum space is all the same 
as our previous work (sec Fig. 8 in Ref. 36). The two di- 
mensional momentum space is thus divided into N x M 
control regions, each with the same energy and angle in- 
tervals. The k-grid point of each control region is chosen 
to be the center of the region: 

\/2m*En{cos0m, sin 0rn. ) , (A3) 



with En = {n + 1/2)Ae and 9m = mAO. n = 
0, 1, • • • , iV - 1 and TO = 0, 1, • • • , M - 1. with the trun- 
cation energy Ecut = Em and 6m = — 1)2tt/M . 

Unlike the electron-LO phonon scattering where the 5- 
function in the scattering is used to carry out the integral 
of k', more specifically ^k', with k' = k — q, here the 5- 
function is used to perform the integral of with 

= J C-^i^r - ■ (A4) 
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